The influence of proximity effects in inhomogeneous electronic states 
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We describe the high-temperature superconductor I^-xSr^CuCU in the underdoped regime in 
terms of a mixed-phase model, which consists of superconducting clusters embedded in the anti- 
ferromagnetic matrix. Our work is motivated by a series of recent angle-resolved photoemission 
experiments, which have significantly enhanced our understanding of the electronic structure in 
single-layer cuprates, and we show that those results can be fully reproduced once a reasonable set 
of parameters is chosen and disorder is properly taken into account. Close to half- filling, many 
nominally hole-rich and superconducting sites have comparatively large excitation gaps due to the 
ubiquitous proximity of the insulating phase. No other competing states or any further assumptions 
are necessary to account for a satisfying description of the underdoped phase, including the notori- 
ous pseudogap, which emerges as a natural consequence of a mixed-phase state. Close to optimal 
doping, the resulting gap distribution mirrors the one found in tunneling experiments. The scenario 
presented is also compelling because the existence of an optimal doping is intrinsically linked to the 
presence of electronic inhomogeneity, as is the transition from the highly mysterious to the rather 
conventional nature of the metallic state above T c in the overdoped regime. The phase diagram 
that emerges bears strong similarity to the canonical cuprate phase diagram, it explains quantita- 
tively the rapid loss of long-range antiferromagnetic order, and also provides for a "superconducting 
fluctuating" regime just above T c as recently proposed in the context of resistance measurements. 
Implications for manganites are briefly discussed. 

PACS numbers: 74.72.-h, 74.20.De, 74.20.Rp, 74.25. Gz 



I. INTRODUCTION 



After almost two decades of intense research, many 
fundamental issues in high-temperature superconductiv- 
ity remain unresolved^. Besides from elucidating the 
mechanism responsible for pairing - although antifer- 
romagnetic (AF) spin fluctuations have emerged as a 
strong favorite^ - there is still no satisfying understand- 
ing of the particularly "tricky" underdoped phase, which 
is most of all characterized by the notorious pseudo- 
gap (PG) in the electronic excitation spectrum. Re- 
cently, remarkable experiments have presented a much 
clearer view of this regime, in particular angle-resolved 
photoemission (ARPES) in La 2 - a; Sr a ;Cu04 (LSCO)£ as 
well as in Ca2- :E Na a ;Cu02Cl2 (Na-CCOC)^ and scanning 
tunneling microscopy (STM) in I^S^CaC^Og+s (Bi- 
22 12 )£&£&£, the latter revealing a granular, inhomoge- 
neous electronic structure on a nanometer scale. Those 
two experimental techniques would neatly complement 
each other - one mapping out the electronic structure in 
momentum space, the other one in real space - were it not 
for the inconvenience that they cannot be performed on 
the same compounds and/or doping levels. Here we want 
to focus on the interpretation of ARPES experiments, 
since the corresponding numerical calculations are not as 
much affected by finite-size effects than those with regard 
to STM data. The ARPES data themselves were shown 
to have provided strong evidence towards the description 
of the underdoped state in terms of a spatially inhomoge- 
neous mixture of hole-rich (i.e., SC) and hole-poor (i.e., 
AF) clustersiS. These conclusions are backed up by a 



variety of other, unrelated, experiments such as Raman 
scattering^ and /iSR-spcctroscopy in the same reference 
compounds (LSCO), as well as by transport measure- 
ments in LSCO and YBa2Cu30jfi£, amongst others. 

Most importantly, the aforementioned photoemission 
(PE) data have granted a unique insight into the elec- 
tronic structure of LSCO over the whole doping ranged, 
and have elegantly highlighted the fundamental differ- 
ences between the underdoped and the overdoped regime. 
Whereas the latter shows fairly conventional behavior - 
comparably to a common, phonon-mediated, supercon- 
ductor - the underdoped phase shows dramatic devia- 
tions from the usual Fermi-liquid physics even in the SC 
state: the spectral intensity was found to have a two-peak 
structure in the anti-nodal direction, across the whole 
range from zero to optimal doping, x c . One of those 
peaks occupies the same energy position (-0.6cV) as in 
the undoped Mott insulator and thus can safely be iden- 
tified with the lower Hubbard band (LHB), whereas the 
second one appears much closer to the Fermi level Ep 
and does so only for finite doping concentrations. In ad- 
dition, finite (although relatively weak) spectral weight 
along the BZ diagonal exists even deep in the insulat- 
ing phase (x=0.03), suggesting that the Mott insulat- 
ing phase contains aspects of metallicity. Very recently, 
vastly improved PE data have become available^. Al- 
though they confirm that the second anti-nodal peak in 
the spectral density distribution is found at (binding en- 
ergy) <x>~-0.2eV, the accompanying distribution is a very 
flat one, and some weight appears close to Ep, precisely 
what should be expected for a sufficiently disordered sys- 
tem. Those measurements also confirmed that the main 



intensity peak itself travels towards E-p as more charge 
carriers are added, and eventually assumes more quasi- 
particle (qp) characteristics, i.e., the associated spectral 
weight distribution becomes sharper. 

Overall, the transition from the undcrdoped to the 
overdopcd regime appears to be a smooth one, and 
seems facilitated by the "fading away" of the LHB 
peak and the concurrent strengthening of the spectral 
intensity close to Ep, both in the nodal and in the anti- 
nodal direction. As was demonstrated before^, most 
of those results can be understood in a rather simple 
framework of competing AF and SC clusters. Doping 
therein amounts to shifting the relative weight of either 
phase; at optimal doping AF clusters will have vanished, 
leaving behind a homogeneous SC phase. It was also 
remarked that the inhomogeneous state in undcrdoped 
LSCO has features of a "ghost" FS in the sense that 
there exists zero-energy spectral weight for anti-nodal 
momenta as well. Those data, once combined with the 
ones from the BZ diagonal, allow for a reconstruction of 
the underlying FS in both theory and experiment. In 
fact, it was shown for LSCO - but not for Na-CCOC - 
that a FS constructed in this way even fulfills Luttinger's 
theorem^. One of the main goals here is to expand on 
our previous calculations by investigating the tempera- 
ture dependence of the mixed-state as well as to explore 
the corresponding phase diagram. So far, LSCO ARPES 
data have only been provided for a fixed temperature 
of T=20K. Higher temperatures are difficult to achieve 
for technical reasons, but nonetheless, particularly in 
the undcrdoped regime, one can make predictions about 
the mixed AF/SC state that should be experimentally 
verifiable. 

A quite puzzling result, however, is the development 
of the doping-induced signal with its observed shift from 
0.2eV below E-p to -0.03eV in the anti-nodal direction. 
The latter value is certainly in line with what one would 
expect for a superconductor with a T c of about 35K as is 
LSCO, whereas the former one stands in stark contrast to 
the small transition temperatures (actually T c =0) in this 
regime of the LSCO phase diagram. This curious aspect 
has not been reproduced or even investigated in depth in 
the previous effort, although some possible explanations 
were suggested to explain the apparent large excitation 
gap in the very underdoped regime. The simplest one is 
that the superconductor enters a strongly coupled regime 
as LSCO becomes more undcrdoped, with the effective 
qp attraction growing. Nevertheless, this is more of a 
restatement of the observed phenomenon, and not an ex- 
planation per sc. Although it seems conceivable that the 
charge carriers lose kinetic energy as the superconductor 
approaches the insulating limit, this viewpoint is diffi- 
cult to reconcile with the presented scenario, because in 
the strong-coupling picture doping leads to a global loss 
in correlation energy, and not to a local change by re- 
placing insulating sites with metallic ones as proposed 
here. The strong-coupling scenario has been amply dis- 



cussed in attempts to explain the PC in terms of pre- 
formed pairs and (classical) phase fluctuations, which 
are of fundamental importance for a strongly coupled 
two-dimensional (2D) superconductor. However, it needs 
to be noted that such thermal phase fluctuations can in 
fact be studied by Monte-Carlo (MC) techniques with- 
out any further approximations^, and the results suggest 
that the increase in the excitation gap is not a strongly 
coupled phenomenon: Since even the very underdoped 
cuprates have spectral weight at Ep along the BZ di- 
agonal, whereas a superconductor with nearest-neighbor 
(n.n.) attraction undergoes a change in symmetry from 
d-wavc to s+id for a sufficiently strong interaction, and 
which in turn is observable by a full gap in the excitation 
spectrum, LSCO apparently does not enter a strongly- 
coupled regime. 

Alternatively, one may assume that at least one other 
order parameter (OP) is present that has not yet been 
accounted for in the previous considerations and which 
would contribute to and enhance the SC gap. One of 
the best known ideas in this this context is the d-density 
wave (DDW)i&. Since it is well known that a n.n. at- 
traction may give rise to instabilities other than d-wave 
superconductivity, it is not unreasonable to assume that 
other orderings may be present as well. The best can- 
didates are charge density waves (CDW), which are well 
known to compete with superconductivity in n.n. mod- 
els. The DDW is a Cooper-pair CDW, and, as proposed 
by Chakravarty et a/., competes with the d-SC phase in 
the same way as the insulating phase. Of course, it also 
has the same momentum dependence as the d-SC, and 
could appear just as such, if experiments not susceptible 
to magnetic effects are performed. Following the spectro- 
scopic data, however, the DDW OP must exist 'on top' 
of the superconductor (i.e., the hole-rich clusters have 
both a finite DDW as well as d-SC OP) and weaken with 
the addition of holes although the local density stays the 
same. There is certainly no compelling reason why this 
ought to be the case. Here, we will propose a much sim- 
pler explanation, with the particular advantage that no 
additional degrees of freedom or exotic phases need to be 
introduced and the desired effect will simply result as a 
consequence of mixing two different electronic orderings 
with vastly different energy scales. 

This paper is structured as follows: In Section II, we 
introduce the Hamiltonians and the mean-field method 
that we use to obtain the results in later sections. Then, 
in Section III we discuss the significance of chemical po- 
tential pinning and contrast that with some calculations 
where disorder effects have been discarded, and fluctu- 
ating stripes emerge as the hallmarks of doped Mott- 
Hubbard (MH) insulators. In Section IV the spectral 
function for mixed AF/SC systems is calculated, and in 
particular we focus on its temperature dependence, to 
complement the previous effort, where only the ground 
state has been studied. In this section, we will also take 
a closer look at the excitation gap in the underdoped 
regime and suggest a simple explanation for the observed 



shift of the low-energy band in the anti-nodal direction 
by showing that this is a consequence of the spatially in- 
homogeneous AF/SC mixture. In the following Section 
(V), we discuss the percolative aspects of the SC/non-SC 
transition and present the phase diagram as calculated 
within the mixed-state scenario and stress its similarity 
with the canonical cupratc phase diagram. The conclu- 
sions, Section VI, summarize our findings and provide an 
outlook for future work, especially with regard to man- 
ganitcs. 



II. THE HAMILTONIAN AND THE 
SELF-CONSISTENT SOLUTION 

Details of the calculations have already been provided 
in the previous work, so here we just repeat the essen- 
tial steps and point out some differences as we consider 
finite temperatures. We investigate the effective mean- 
field Hamiltonian i?HF, 
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that can be derived from the original extended Hubbard 
model in a Hartree-Fock approximation, assuming both 
finite SC and magnetic order exist in certain parts of 
the sample, c- in Eq.Q) are electron (creation) opera- 
tors on a two-dimensional quadratic lattice of N=LxL 
sites, and t is their concomitant hopping amplitude be- 
tween n.n. sites i, j. We chose t=l as the energy unit. 
The particle density ni=^2 
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is determined by the 
chemical potential \i\ and both are assumed to vary spa- 
tially. In Eq.Q), we have also introduced the local spin 
operator s\ =i;(n\i l -n\i) . The local magnetic OP m\ and 
the site-dependent (d-wave) SC OP Ay are related to the 
fermion operators, 
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where (. . .) denotes expectation values after thermal av- 
eraging. Vy (defined on links i,j) is the n.n. at- 
traction that appears in the original extended Hub- 
bard Hamiltonian and it is responsible for the occur- 
rence of SC (or charge-order (CO), if a charge-density 
wave were to be considered), whereas Ui is the on-site 
Coulomb interaction, leading to AF order. The Hamil- 
tonian Hhf=H^ if +H^ f is quadratic in electron opera- 
tors, which makes it a single-particle Hamiltonian with 
27V basis states, that can be readily diagonalized using 
library subroutines. H^ F (third line in Eq.Q) is a solely 
classical energy term, with no operators involved. It sim- 
ply changes the configuration energy, but will not be of 



importance for the rest of the paper. Eqs. Q), (J2J) con- 
stitute a self-consistent system, which can be solved it- 
eratively. This amounts to performing a Bogoliubov-de 
Gennes (BdG) transformation22i2i24, whereby the orig- 
inal operators c\ a are expressed in terms of new Bogoli- 
ubov quasiparticlcs 7„ CT : 

N 
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This, together with Eq. J2J, leads to the following ex- 
pressions for the OP's as functions of the wave functions 
a„(i), 6„(i): 
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where (7i cr 7 nCT )={l + e /3En, '}~ 1 , i.e. the Fermi function 
/ with properties f(x)=l-f(-x), (3=1 /T the inverse tem- 
perature and E na are the eigenvalues of Eq.Q). Whereas 
previously we had focused on ground state properties 
only, this time we will investigate specifically the tem- 
perature dependence. At the start of the iteration initial 
values A?= , m? arc chosen, and the resulting Hamiltonian 
Eq. (Q is diagonalized, leading to wave functions a°(i), 
&°(i) and eigenvalues E® lcr . Those are then used to com- 
pute A«, m\ via Eq.QJ, and so forth. Those iterations 
are stopped once ^^Znl Ajj-Ajj x |<e, where e typically 
is 10~ 5 , and a similar condition applies to the magnetic 
OP as well. In case no magnetic order is present a n =b n , 
leading us back to the original BdG equations. 

Using mean-field approximations, it is very well 
known that for Hamiltonians such as i?HF chargc- 
inhomogeneous regions appear spontaneously in the low- 
doping limit (<5=1--^ J7 n i> <^1) m the form of stripes^. 
This regime is difficult to handle numerically, and may 
not even fully reflect the physics of underdoped cuprates 
since it neglects the influence of chemical disorder. Natu- 
rally, one would expect such terms due to the presence of 
the Mott insulating phase and its associated poor screen- 
ing properties. Therefore, in our case doping is always 
accompanied by a change in the local chemical potential 
/iti, which is set so that the resulting local charge den- 
sity rii is the one of optimal doping. This is achieved by 
introducing "plaquettes" , which are small patches of re- 
duced chemical potential (and thus charge) and enhanced 
pairing^. This expresses the view that Coulomb po- 
tentials, maybe in conjunction with the AF background, 



partially bind holes to their donor. Then, the chemical 
potential is composed of a local term and a global one, 
Ho; the latter still has to be determined self-consistently 
to obtain the desired global charge (or hole) density. In 
order for our results to reflect the PE data, we will com- 
pare the computed values of /xo with the ones obtained 
via PE in LSCO. 

From the knowledge of the eigenfunctions, all relevant 
observables such as the single-particle spectral function 

A(k,u>) = — JmG(k,w), (5) 



G(k, «)=£,,,_ 



,ik(i- 



function, and the local density of states 
N(i,uj) = —ImG(i,u> 



c icr c i'cr) being the usual Green's 
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can be calculated in a straightforward fashion. The den- 
sity of states N(ui) = jj^2 i N(i,oj) follows immediately. 
Those observables are the crucial quantities with regards 
to ARPES and STM experiments, respectively. In the 
two limiting clean cases, {n}=l and with no SC regions 
present, as well as (n)=0.75 (=nsc) and without resid- 
ual AF islands, the BdG equations converge fairly quickly 
and one obtains the curves A(T) and m(T) (spatial in- 
dices are not necessary since the systems are homoge- 
neous), which can be compared with results from solving 
the self-consistent equations in momentum-space. The 
temperatures T c upon where A(T c )=0 (m(T c )=0) de- 
fine the corresponding critical temperatures. For the 
two cases considered, one finds T^ F ~1 and T c sc ~0.35 
((n)=0.75) for the parameter values U=5 and Vy=V=l, 
respectively, which we will adopt from here on, unless 
mentioned otherwise. It is mostly important for our pur- 
poses that they are sufficiently separated, but we do point 
out here that the separation of energy scales is not as 
large as in LSCO, where T^ F RilOT^ c , an issue we will 
return to later. 



III. THE CHEMICAL POTENTIAL AND THE 
INFLUENCE OF DISORDER 

As mentioned above, the local chemical potential is 
chosen so that the resulting charge density at hole-rich 
sites stays approximately constant, and further doping 
is mimicked by introducing more such charge-depleted 
regions. Effectively, as we will also show below, this re- 
sults in a pinning of the global chemical potential [Aq at 
the original Fermi level. Since this pinning, according 
to ARPES, holds up to roughly optimal doping upon 
where a steep drop in fio in agreement with common 
band-structure theory is registered, it is natural to as- 
sume that the chosen charge density level for hole-rich 
clusters should be the one of optimal doping and we ad- 
just fa accordingly. 

FigQ] shows the resulting low-temperature N(cu) as a 
function of density, obtained after averaging over 3 con- 
figurations on a 40x40 lattice. Given the comparatively 
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FIG. 1: Density of states N(u) for the mixed SC/AF sys- 
tem on a 40x40 lattice at low temperatures. Parameters are 
given in the text. The low-energy peak at small x is not intrin- 
sically connected to the prominent coherence peak observed 
for 2=0.20, but is an accidental consequence of the plaquette 
structure chosen. The inset shows the chemical potential, /io, 
as a function of the overall density, no is almost completely 
pinned to the original Fermi level uj—0 until the mixed-phase 
regime is replaced by a homogeneous (SC) state, upon where 
fio drops in accordance with band-theory. 



large size of our lattice, even such a relatively small 
number should be sufficient to guarantee configuration- 
independent results. By comparing results obtained for 
individual configurations we also have confirmed that 
there are only minor variations in between them and the 
systems therefore are close to self-averaging. The large 
spike associated with the MH insulator appears at wmh~- 
2, and is determined by our choice of U=5. Upon doping, 
in-gap (i.e., between this MH gap) states are created, and 
they are continuously filled as the system evolves towards 
x c . The MH peak in N(tu) diminishes in size, and is even- 
tually replaced by the SC coherence peak positioned at 
w~-0.7. Nevertheless, as the original gap is filled by the 
hole-rich regions, the spike at wmh remains in its original 
position throughout the doping process. This is the LHB 
pinning observed in ARPES. The small signal at wmh/3 
((?i)=0.95) corresponds to the second band found along 
the anti-nodal direction upon doping (see Refs. I3l fl5f> . 
Whereas in PE, however, this peak moves towards Ep, 
in FigQ]the opposite behavior is observed. This issue has 
caused some confusion in the previous worki£, and it will 
be discussed in great detail in a later section. 

The inset in Fig ^ shows how ^o develops: it remains 
pinned at the original level c<j=0 for densities down to 
ngc; but drops steeply once the system has become ho- 
mogeneous. Although this is congruent with what has 
been found in LSCO, by no means does this describe all 
underdoped cuprates as a more conventional behavior is 
found for Na-CCOC, where the (surprisingly) broad band 
associated with the AF insulator moves continuously to- 
wards Ep with the addition of holes. It is also important 
to consider that the scenario presented here could still 
be somewhat oversimplified, since a careful analysis of 
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FIG. 2: N(w) derived from a MC integration of Fhf on a 
16x16 lattice. J~Um— 0.6 and T=0.04. Note that these are 
original data without Lorentzian broadening. The inset shows 
the particle density as a function of the chemical potential. 



spectral data seems to suggest that although an interme- 
diate band is indeed created upon doping, the chemical 
potential does move slightly downward inside this band 
as more holes are added. Whether the observed small 
shift in [i between x=0 and a;=0.25 (see inset FigQJ de- 
scribes the same effect is, as of now, not clear, but we 
will further comment on this below. 

The importance of disorder and its correlations with 
the electronic structure were recently demonstrated in 
STM investigations of slightly undcrdoped Bi-2212±&, 
and the influence of dopant atoms was related to a change 
of the local chemical potential as far as the CuO-layers 
are concerned. There, it was also remarked that the 
charge fluctuations are surprisingly small, of the order 
of 10%. This is basically the viewpoint taken here as 
well. Alternatively, the case where doping is not ac- 
companied by disorder should also be considered, since 
there are indications that not all high-T c compounds are 
equally disordered. Although in principle doping with- 
out accompanying change in the local chemical potential 
lii (as described in Section III) can be achieved by sim- 
ply reducing /xo, in practice this is not possible within 
the mean- field framework employed here. As is well 
known, Eq.Q) is unstable towards phase-separation in 
the form of stripes, and the intermediate regime between 
the Mott insulator at (n)=l and the stripe phase at a 
much smaller, precisely determined, value of (n) cannot 
be accessed (in particular for the lattice sizes considered) . 
Yet, as it turns out, this phase can be adequately de- 
scribed via a MC technique that was introduced in the 
previous paper—, albeit at the price of being limited to 
smaller lattices. N(uS) resulting from a MC integration 
of Hhf at low, but finite, temperature is displayed in 
Fig|3for a relatively large lattice of 7V=16xl6 sites. For 
this calculation spatial variations of ^i and Ui were not 
considered. Doping with charge carriers is mimicked by 
reducing the (global) chemical potential, as underlined 
in the inset of Fig(2j As [i is reduced, the density drops 
concurrently. This process, however, may not be entirely 



smooth: in fact, a tendency to phase separation is found 
around /x=-0.5, between densities (ra)=0.94 and (n)~0.87 
(x~l/8), signaling a first order phase transition. Similar 
results were also found for models that implement the 
SO (3) symmetries of the classical spinal, which shows 
that the Ising spin approximation employed here is not 
a crucial one. Even in the slightly doped case, the MC- 
averaged density is constant throughout the lattice, but 
snapshots show that single sites are either almost per- 
fectly half-filled or, to a much smaller degree, possess a 
density (n}~0.63, this value mainly being a function of 
the magnetic coupling J=^Um, which, instead of U, en- 
ters the MC process^i. Hole-rich sites sometimes arrange 
themselves along straight or diagonal lines, without form- 
ing full stripes. Essentially, the MC process describes a 
phase of dynamical stripes, which are broken in pieces ei- 
ther because of thermal fluctuations or because the den- 
sities involved do not allow for a full stripe to be formed. 
Notably, in this regime the original spike in N(ui) at en- 
ergy J remains at its position (FigEI long-dashed line) 
and is only diminished in size, whereas at the same time 
finite spectral weight appears at low-energies, w~J/3- 
J/2. One peak clearly relates to the fully occupied sites, 
whereas the latter one reflects the residual magnetic gap 
at hole-rich sites, a view underscored by a direct mea- 
surement and also by the observation that its position 
scales with J over a range of J (0.3< J<0.8) values ob- 
served. This is quite similar to what is found in Fig^ but 
upon further doping this two-peak structure is lost, and 
replaced by a smooth function without any noticable fea- 
tures in N(u>), even though charge-disorder prevails. For 
example, even at the stripe-favoring hole concentration 
level x~l/8, do we not find the naively expected stripe 
state, but instead stripe-like states with charge distribu- 
tions centered around densities ni~0.94, ri2~0.75, and 
similar results for other values of \i. In other words, we 
have disorder without pinning. This stands in contrast to 
LSCO PE data, where pinning and finite excitation gaps 
arc observed across the underdoped regime. Similar re- 
sults were found for a wide variety of values of J, and thus 
they should be regarded as representative for disorder- 
free Hamiltonians like Hhf and not as accidents caused 
by unfortunately chosen parameters. All this suggests 
that a certain amount of disorder has to be present to 
account for the observed stationarity of the MH-band in 
LSCO, whereas for Na-CCOC possibly the weak-disorder 
scenario is the more appropriate one. However, because 
an inclination to chemical potential pinning does exist 
(Fig(2Jl, disorder strength may not be overly large to in- 
duce the effects described here. Note, however, that at 
this point it is difficult to decide whether such solutions 
simply result from the original mean-field approximation 
or represent genuine low-energy states for doped MH in- 
sulators and further investigations may be needed. 



IV. SPECTRAL FUNCTIONS IN THE 
PRESENCE OF COMPETING STATES 



The BdG equations were solved on lattices of the size 
iV=40x40, which allows for a detailed comparison with 
experimental data. Although larger systems could in 
principle be used, the great number of iterations nec- 
essary to achieve self-consistency at temperatures not- 
too-far from clean-case critical temperatures renders such 
systems practically useless. We have also observed that 
as the lattice is increased the convergence at high tem- 
peratures suffers and additional iterations are required 
to arrive at a self-consistent solution. As in the previous 
effort, charge-depleted "plaqucttcs" were distributed on 
the lattice in a random fashion. Those plaqucttcs were 
allowed to overlap, but care was taken to ensure that 
the area covered (i.e. the relative amount of sites with 
reduced chemical potential) grew linearly with doping. 
Since the overlap probability increases with the number 
of plaquettes already present, their number has to grow 
supralincarily with the hole concentration. In this fash- 
ion, however, truly random configurations can be pro- 
duced, and the corresponding self-consistent equations 
solved. 

One such example can be seen in FigQJa)-(c), where Ai 
and m,i are shown for hole concentration a;=0.15 at three 
different temperatures, starting from T~0 (FigQJa), 
(af )) to a finite temperature T<T^ C ((b), (bl)) (see Sec- 
tion V), and finally at TwT c sc ((c), (cl)). Given out- 
values for ngc=0.75 and ?iaf=1, about 40% of the to- 
tal area arc insulating. Dark color translates to large 
OP amplitude, chosen in a way so that the darkest pos- 
sible ones coincide with the corresponding amplitude in 
the clean case. As far as the SC OP is concerned, an 
increase in temperature at first transforms sites at the 
boundaries into metallic ones, but leaves the OP strength 
inside larger clusters unaffected (Fig^Jb)). Eventually, 
this leads to a breakup of the original SC "backbone" 
into a set of single, isolated islands (FigQJc)). Those 
small clusters, however, arc capable of surviving to sur- 
prisingly large temperatures, in this particular case up 
to T-0.28 (3/4T c sc ). The local AF order, on the other 
hand, appears more stable as it persists almost up to T^ 
with only a small loss of overall volume (Fig0Jcf )), and 
thus there exists an extended temperature range where 
metallic and AF regions coexist. This observation can 
be made regardless of the overall amount of insulating 
clusters. The phase diagram that can be derived from 
data like those will be unveiled in Section V. 

The choice of plaqucttcs, i.e., their size, does have cer- 
tain effects on the outcome, and this is demonstrated in 
Figs|3fa), (b), which shows m\ for (a) the procedure used 
throughout the paper, and (b) where a large number of 
very small (i.e., of size 1) "plaquettes" were used— In 
the latter case, the clusters appear much more irregular 
and fractal, as expected. Although the resulting spectral 
functions are hard to distinguish, and thus our results 
not obviously affected by any particular choice, it does 







FIG. 3: Influence of the plaquette geometry on the local AF 
OP distribution, (a) corresponds to the procedure of using 
extended plaquettes, mainly adopted in this work, whereas in 
(b) the size of the plaquettes is 1, leading to a more "fractal" 
shape of the resulting clusters. Particle density (n)=0.85 and 
T~0 in both cases. 



influence the OP distribution function as we will discuss 
in a later section. 

The spectral functions related to mixed AF/SC states 
such as in FigQ] are shown in FigJSJ those for the two 
"clean" cases have already been presented in the previous 
paper— and are therefore not displayed here. There, it 
was also demonstrated that the spectral functions in the 
mixed state can effectively be obtained by superimposing 
the spectral functions of the respective clean phases on 
top of each other, provided the two signals are weighted 
by the relative amount of either phase. Thus, at the low- 
est temperatures considered, both the two-peak structure 
at (tt,0) and (although it can be barely seen here) at (tt/2, 
7r/2) is visible as is the presence of a small gap at E-p in 
the anti- nodal direction (FigJSfa)). The spectra them- 
selves are "smeared out" and poorly defined, especially 
in the anti-nodal direction, similar to what is observed 
in ARPES. The AF branch is found at w~-2, consistent 
with Fig^ whereas the maximum of the low-energy in- 
tensity is at w~-0.5, but the corresponding distribution 
is very broad and intensity is found almost up to i?p— 
As T is increased and SC states become metallic ones 
(see Fig0J, one observes the accumulation of spectral 
weight right at E-p (FiglSJc)). This takes place in the 
temperature window somewhat (but not too much) be- 
low T C 5C , and the FS crossing becomes much better de- 
fined for higher temperatures (Fig^d)). Nevertheless, 
remnants of the insulating regime still survive at those 
temperatures, and the two-peak structure in A(k, lu) is 
preserved up to almost T^ F . Note, however, that the 
SC phase is replaced by a metallic one and not by the 
competing AF state. The latter scenario is certainly also 
possible, particularly in the clean case, and if a first-order 
phase transition is separating both phases. In such a case, 
the signal close to E-p would actually vanish at elevated 
temperatures, and the only remaining gap would be the 
one corresponding to the LHB. Unfortunately, ARPES 
as of now, cannot tell us what is happening since mea- 
surements in this temperature range are currently not 
possible^. As long as disorder is important, however, we 




FIG. 4: Order parameters Ai (left column) and m\ (right 
column) at different temperatures T~0 ((a), (al)), T=0.10 
((b), (bf)), and T=0.25 ((c), (cl)) at overall density (n)=0.85. 



think that the presented scenario is the correct one. The 
position of the LHB in any case is strongly influenced by 
the temperature as well, quite unlike what is observed if 
the carrier-density is changed. It shifts to lower binding 
energy, and weakens at the same time. 

At this point we want to emphasize that the transi- 
tion between different regimes is a very "soft" or gradual 
one, and particularly the loss of long-range SC order is 
not coinciding with any dramatic changes in the spectral 
functions, unlike in the case of the common BCS transi- 
tion, which entails the sudden emergence of Bogoliubov 
quasiparticles. Fig \5\ testifies to that. With the demise 
of superconductivity one should, of course, expect the 
emergence of a complete FS somewhere below T c s c , but 
this is not what is occurring here, at least in the weakly 
doped regime: since only a small volume percentage is af- 
fected by the SC/metal transition, and with most of the 
spectral weight residing in the AF branch, the changes 
in the spectral functions are minuscule and they do not 
resemble those of a true qp. Clearly, the properties of 
the sample are dominated by the presence of insulating 
regions. All this is quite typical for systems close to a 
percolation threshold, where a small change in volume 
fractions (which arc effectively measured by ARPES) can 
induce a drastic change in transport properties. 
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FIG. 5: A(k, lo) at density (n)=0.85 at 4 different temper- 
atures: T~0 (a), T=0.I0 (b), T=0.30 (c), and (d) T=0.80. 
Even though a FS is established for k~(7r,0) at intermedi- 
ate temperatures, the SDW branch survives up to the largest 
temperatures and is responsible for the depletion of the low- 
energy density of states. 



A. Proximity effects in inhomogeneous states 

There has been considerable progress in the quality 
of data provided by PE experiments since the first re- 
sults have become available. Whereas originally ARPES 
results only revealed the large excitation gap of about 
we- 0.2 eV (a;=0.03) in the anti-nodal direction, later re- 
sults have shown that spectral weight can in fact be found 
for all energy values between wg and E$^. Nevertheless, 
it was confirmed that the main intensity indeed resides at 
we; this main intensity peak moves towards E-p at subse- 
quent doping, but remains poorly defined, and becomes 
qp-like only close to optimal doping. Fig^does not fully 
reproduce this behavior, since it suggests that the gap is 
actually increasing with doping rather than decreasing. 
Certainly the initial large value of u>e can be described 
in terms of a very strong attraction or an additional, but 
hitherto hidden, OP such as DDW, by making V a func- 
tion of x; nevertheless, in the scenario where doping is 
achieved by replacing hole-poor sites with hole-rich ones, 
it is mysterious as to why those parameters are changing 
even though the (local) density remains constant. Here 
we attempt a different explanation for the shrinking of 
the gap, and one that does not rely on additional assump- 
tions beyond those of the existence of a mixed-state. 

To do so we first remark that the choice of parameters 
in Ref[l5| (as well as in this work) has been somewhat 
unfortunate as the distinction between the AF and the 
SC energy scales that is present in the cuprates is lost to 
some degree, for V—U/b is unnaturally large, at least as 
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FIG. 6: The spectral functions for a mixed AF-metallic model, 
with U—5 and l/=0.0 and densities (a) (n)=0.95, and (b) 
(n)=0.90. Even though the metallic clusters are ungapped 
courtesy of Hhf, there is yet a finite excitation gap in the 
anti-nodal direction close to half-filling ((a)), caused by prox- 
imity effects. The same data as in (a), but for V=0.6, are 
plotted in (c) and (d) is the spectral function for the homoge- 
neous superconductor with V=0.6 for (n)=0.75. The shrink- 
ing of the gap close to (n, 0) between (c), (d) is obvious. 



far as LSCO is concerned. The self-consistent solution 
of Eq.Q tells us that in a mixed-phase sample due to 
proximity effects there are always sites or regions that 
carry induced OP amplitudes, namely those at or close 
to phase boundaries (Fig^J. Therefore, sites that are 
supposed to be solely SC, due to the local term in Hhf, 
may have a finite magnetization as well and vice versa. 
To better understand the relevance of such effects, Hhf 
was solved again, but this time with V=0 in order to 
disentangle the AF from the SC OP. This is important 
because in actual materials such proximity effects, due to 
the large energy associated with AF ordering, may over- 
shadow the influence of the SC term and actually deter- 
mine the observed gap. In contrast, for the particular 
choice of V/U made here as well as in the previous work, 
the "clean" SC gap Agi=y • Ai is about as large as the 
magnetic gap A m =^-mi on sites i located next to an AF 
cluster. Therefore, we reduce V in an effort to correctly 
describe the properties of undcrdoped cuprates, starting 
out with the simplest case, V=0. The spectral functions 
now emerging arc shown in FigsEfa), (b). Even though 
the low-energy intensity is not very significant, one can 
conclude that most of the low-energy weight in the anti- 
nodal direction is concentrated at binding energy uj^-0A 
for x=0.05 and appears to have a finite gap, although 
the corresponding interaction V has been set to zero! An 
inspection shows that this finite gap is caused by prox- 
imity effects, which dominate for small doping levels. As 
doping progresses, this intensity travels further towards 



Ep, and actually starts to cross Ep for 0.05<a;<0.10. Al- 
though the signals are blurred due to the inhomogeneous 
nature of the sample, the main issue, the existence (or 
enhancement) of excitation gaps close to (n)=l is clearly 
demonstrated in Fig|^a). We also want to reiterate here 
that the same "blurriness" in the spectral intensities has 
been observed for samples with intermediate (a;=0.07) 
hole concentration^. 

Furthermore, similar conclusions can be achieved by 
directly measuring the AF OP distribution, N(m), which 
is depicted in Fig[7| for £=0.05, and which serves to 
substantiate and clarify our arguments. For such a 
small number of holes, N(m) has the expected peak at 
w m ~0.75 (thus, A m «2, see Fig^l in parallel with the 
one found for half-filled systems and stemming from sites 
that are located deep in AF territory, plus low-energy 
contributions below cl>«0.3. They stem from U=0 sites 
in Eq.Q, that suffer from proximity effects. For x=0.05 
there are only few charge-depleted, uncorrelated clusters, 
and many such sites are in close proximity to hole-poor 
ones. Therefore, the low-energy distribution is peaked 
at a comparatively large value, u>i c ~0.12, although the 
distribution itself is somewhat broad (see FigEJa)). This 
distribution, however, depends to quite some degree on 
the specific shape of plaquettes chosen, which is demon- 
strated in Figs[7fb), (c). The latter one, corresponding 
to a OP distribution such as depicted in Fig|2{b) shows 
that for a sufficiently fractal configuration, the AF OP 
can assume any value between zero and w m , whereas the 
former one resides between the two extremal cases con- 
sidered here, and has a distribution peak at wi o ~0.2. Al- 
though this ratio is closer to the actual ratio 1/3 (0.2cV 
vs. 0.6eV), our arguments are not affected by any par- 
ticular choice, and we will mostly stick to our original 
geometry. 

In any case, adding holes changes N(m) considerably, 
as it becomes dominated by non-interacting sites and the 
peak of the low-energy part of N(m) is shifted towards 
uj\ c =Q. This is demonstrated in Fig|Hl which shows the 
N(m) at different densities as one progresses towards 
x c . Whereas initially (FiglHJa), same as FigEJa)) the 
low-energy part is dominated by finite energy contribu- 
tions, subsequently the main peak in the distributions 
shifts to smaller energies (Fig|SJb)). Thus, at the onset 
of doping, the majority of states has a larger excitation 
gap than predicted by the value of V (which is actually 
zero here), and it is those states that have the largest 
weight and, therefore, intensity in PE; still, there are a 
few states with a very small (or zero) mi, and they are 
responsible for the finite spectral weight at Ep (along 
the BZ diagonal). As doping progresses, less and less 
metallic sites will be influenced by neighboring insulating 
ones, and eventually the vast majority of sites will have 
only a very small or even zero magnetization (Fig|SIc)), 
causing the spectral intensity peak to move towards Ep. 
A simple calculation confirms this picture: the average 
magnetic gap A m (counting only the low- energy contri- 
butions) changes from A m =0.33 ((n}=0.95) to A m =0.26 
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FIG. 7: N(m) as a function of energy for the density (n}=0.95 
(i.e., 20% non-AF sites). The main peak at cj m ~0.75 pertains 
to the homogeneous AF phase, whereas the low-energy values 
stem from boundary sites, (a) is the distribution for the 12- 
site plaquette as mostly used in this work, (b) is the 6-site 
plaquette, and (c) the 1-site plaquette. The samples become 
more fractal when going from (a) to (c), and this is reflected 
in the very broad OP distribution function. 




FIG. 8: N(m) at low temperatures as a function of energy 
for varying densities: (a) (n)=0.95, (b) (n}=0.90, and (c) 
(n)=0.85 in the mixed-phase regime. V=0 for non-AF sites. 



(_(n)=0.90), to A m =0.20 ((n)=0.85), and finally to 
A m =0.11 ((n)=0.80). These results refer to FigsEIa)- 
(c) but similar ones were obtained for other configura- 
tions as well and demonstrate that any low-energy gap in 
mixed insulating/metallic systems decreases continuously 
upon replacing insulating with metallic sites. An iden- 
tical behavior was observed in STM as welU&. Choos- 
ing a finite V will not change these results aside from 
adding a small additional gap (if V is small compared 
to U) and preventing A(k, u>) from actually develop- 
ing a FS crossing. For any given site, the local effec- 
tive gap (at least in the mean-field approximation) will 
then be A e ff 



of both AF and SC gaps 



A^, i-e 



it is a combination 
The relative distribution of 



FIG. 9: Relative distribution of the total gap A e g for the 
density (n)=0.80 (i.e. 80% SC) for V=l (a), V=0.80 (b), 
V=0.60 (c), and V=0.60, but a more fractal configuration 
(d). The best agreement with STM data is found for (c). 



A off ((n)=0.80), N(A cS ), is shown in Fig for differ- 
ent values of V. The results for V=l, the value chosen 
in ReflTa, can be found in Fig[5Ja), and most impor- 
tantly, iV(A c ff) is peaked at the comparatively large en- 
ergy w~0.75 (compare, e.g., with FigQJ, and thus As 
dwarfs residual AF gaps (A m ~0.1). In this case, the 
average gap A e g actually increases with hole-doping, as 
relatively small proximity-induced AF gaps are replaced 
by large SC ones. Then, a much better choice to under- 
stand PE results in LSCO is found in V=0.6 (FigG^c)), 
where both Agi, A m i (in hole-rich areas) are of the same 
order, and N(A c g) is peaked slightly below ui=0.2. In 
this case, A e ff decreases with the addition of holes, from 
the original value (x=0.05) A e ff~0.33 (superconductivity 
is almost negligible there) to, eventually, A e ff«0.2, where 
the gap is mostly defined by V. In other words, we argue 
that the observed large effective gap at small hole concen- 
trations is caused by intermixing antiferromagnetism via 
proximity effects, which is then gradually replaced by the 
pure SC gap as doping progresses. This provides in our 
view the simplest and most natural explanation for the 
doping dependence of the LSCO low-energy gap. This 
picture agrees also with the recent notion that this gap is 
composed of a low-energy contribution following d-wave 
behavior, and a larger one, leading to an overall devia- 
tion from the simple cos k x -cos k y formSi. Figs Etc), (d), 
which show the anti-nodal gap for the same configura- 
tion as before, but with V~=0.6, changes from we~0.5 
to we^0.2 between (n)=0.95 and (n)=0.75 confirms this 
once more as does N(uj), where the small original peak 
(Fig|TJ grows and evolves towards E-p (i.e., w~0.2) with 
doping. Yet, we stress that proximity effects are impor- 
tant since without them the peak associated with super- 
conductivity simply would remain in its position. The 
critical temperature for T^=0.6 is T c ~0.14, and thus the 
respective clean case AF/SC critical temperature ratio 
~l/7, reasonably close to the actual LSCO ratio 1/10. 
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Although this leads us to believe that an even more faith- 
ful representation of ARPES results can be achieved by 
choosing a smaller value of V, the resulting difficulties 
in the self-consistent scheme prevent us from doing so. 
Already, however, iV(A ff) in FigHJc) strongly resemble 
the one found via low-bias STM in Bi-22132£: its shape 
is almost Gaussian, with only a minor tail, quite unlike 
Fig|S[a), (b), (d). We also mention here that in RefyfJ 
the gap distribution is mirrored by a similar charge distri- 
bution, as is the case in our model. A logical next step is 
the calculation of iV(i, u>) and a detailed comparison with 
STM data. Unfortunately, this calculation suffers from 
serious finite-size effects^!, and important questions such 
as whether sites associated with very small (large) mi are 
associated with a comparatively prominent (small) co- 
herence peak cannot satisfactorily be resolved. For this 
reason, we forgo a discussion of N(i, u>) for the time be- 
ing. 

Our goal here was to demonstrate that it is necessary 
to carefully chose the parameters to fully reproduce the 
LSCO PE data. Proximity effects ought to be ubiqui- 
tous in mixed-state phases, and it has now been shown 
that their effects have measurable consequences. Neutron 
diffraction experiments by Lake et alm^ showed that the 
AF response in undcrdoped LSCO increases upon the 
application of a magnetic field, which supresscs super- 
conductivity, a result which may be consistent with the 
existence of SC/AF coexisting regions in inhomogeneous 
materials, where the same behavior can be expected. 



B. The Fermi surface 

The corresponding FS can be determined by measur- 
ing the spectral weight at Ep (integrated over a small en- 
ergy window, Ao>~0.2-0.3), and it is shown in Figs llOf al- 
(d) as it evolves with temperature ((n)= 0.90). Close to 
T~0, significant weight (and a qp-like peak) is found only 
in the neighborhood close to the BZ diagonal, as already 
demonstrated in the previous paper—. The resulting FS 
"arc" widens as the temperature is raised, indicating en- 
hanced metallicity, until it becomes a full (electron- like) 
FS at about T=0.80 fFiglWdl). Nevertheless, even for 
the temperatures distinctively above Tj? c , does the FS 
not fully resemble the one of a conventional metal. In 
fact, even for T>2T C SC does the FS appear "blurred" at 
momenta (0, n), (n,0), and only above ~T c Ai? is a reg- 
ular, sharply-defined FS recovered. Thus, if the PG is 
caused by a remnant insulating phase, it must at the 
same time be characterized by a poorly defined FS at 
the anti-nodal points. Similar results arc found for other 
densities as well. In particular, this has important impli- 
cations for the SC transition for x>0.10 (see Fig ll2|l ; in 
a conventional superconductor, this transition is accom- 
panied by the restoration of the underlying FS in the 
anti-nodal direction, something that docs not happen in 
the case regarded here. Therefore, the change from the 
SC phase into a "metallic" one will appear as much more 
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FIG. 10: Spectral weight at the Fermi level at different tem- 
peratures for Hhf: T~0 (a), T=0.20 (b), T=0.50 (c), and 
T=0.80 (d). The original "arc" (indicated by white lines) 
expands as the temperature is increased, and becomes a full 
FS for temperatures much larger than T' c f . Even then, be- 
cause of residual insulating clusters, the FS in the anti-nodal 
directions is not particularly well defined. 



gradual, and less abrupt than a conventional SC/metall 
transition. 

In the mixed regime, the length of the FS arc is to 
a large degree independent of the hole density, whereas 
the spectral intensity at Ep grows continuously. This 
observation, made in PE experiments and later repro- 
duced in calculations based on the AF/SC mixed-state 
pictur e) 3 ! 15 , is naturally interpreted as a consequence of 
the increased relative weight of the SC phase, and sub- 
sequently, of the importance of zero-energy excitations, 
within the sample upon doping. Those observations are 
somewhat reversed when the temperature dependence of 
the FS is considered: as the arc is spreading out with en- 
hanced T (T<0.8 for the parameters considered here), 
the intensity (along the BZ diagonal) stays the same 
(within 15%); this so happens because SC regions are 
turned into metallic ones, and the states most affected 
by that are the ones that are close to the BZ diagonal, 
where A is smaller. States along the BZ diagonal are 
not affected since they are already metallic. Neverthe- 
less, the spectral intensity along the anti-nodal direction, 
which is independent of doping to a large degree, shows 
a distinct increase with temperature within the mixed- 
state scenario as some marginally AF sites are turned 
into metallic ones. Note, that it is only important in this 
context that the SC/AF energy scale are sufficiently sep- 
arated so there exists a broad temperature range without 
SC order, but an intact AF phase. This certainly applies 
to the cuprates, and therefore, these predictions are ver- 
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ifiable within PE experiments as long as the necessary 
temperatures can be achieved. 

Furthermore, the scenario presented here naturally ex- 
plains the behavior of the (integrated) spectral weight at 
E-p. It was shown in Refy (Fig. 4 therein) that this quan- 
tity (as well as the qp-density derived from Hall effect 
measurements) grows continuously with doping, which 
naturally reflects the increased amount of metallic areas 
in the sample, as already remarked in Reflla. Yet, the 
rate of increase is comparatively small for a-<0.1, but in- 
creases strongly thereafter, whereas this rate is constant 
following Reflla. Proximity effects provide us with an ex- 
planation: close to the parent insulator, many nominally 
"free" sites are still to some degree AF ordered, and the 
metallic volume is smaller than one would expect accord- 



TABLE I: Integrated spectral intensity near the Fermi level 
as a function of the electron density, (n) . 
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ing to the doping level. To make up for the "lost" sites - 
which has to happen eventually at optimal doping - their 
number has to increase at a faster rate at higher doping 
levels. To underline this point we have integrated the 
near-iJp spectral weight (within a window of Aw=0.3eV), 
again for the case of V=0, and have tabled those values, 
ii_3, in Table I. I\ (I2) refer to different k-space win- 
dows (along the BZ diagonal), whereas I3 refers to data 
obtained by using Aw=0.4, but the same BZ cut as in I\. 
Irrespective of how the window is defined, we find that 
the rate of increase is small for x < 0.05, but considerably 
larger afterwards. It is quite impressive that even such 
details of PE data can be reliably reproduced within the 
simple mixed-state picture. 

Still, even at low temperatures, the high-intensity spec- 
tral weight together with the low-intensity signals in the 
anti- nodal direction fFig llOf a) - ) evokes memories of a 
complete FS, as one would expect it for a non- interacting 
system of ~nsc electrons. Experimentally, it was demon- 
strated that at least in the case of LSCO a FS con- 
structed in this way approximately fulfillcs Luttinger's 
theorem^, and a similar result was found theoretically 
as welli&. This interesting result suggests that the under- 
doped phase is actually closely linked with the weakly- 
interacting metal, and is, for example, not separated from 
it by any FS transition. Our results here imply that the 
FS that ought to emerge as the temperature is raised is 
the one of the underlying, hole-rich metallic phase, which 
in itself is independent of the overall doping level, that is, 
it should be the same FS independent of doping. ARPES 
data at elevated temperatures could confirm this picture. 



V. PHASE DIAGRAM 

Based on the considerations presented above it is 
tempting to develop a generic phase diagram of -Hhf in 
the T-x-plane, across the whole doping range. This is rel- 
atively straightforward in the homogeneous case, where 
the condition A(T c )=0 determines the critical lines. In 
the inhomogeneous scenario presented here, a more care- 
ful approach needs to be taken in order to distinguish 
between global and local phenomena, both of which char- 
acterize the sample as discussed above. In the mean-field 
scheme employed here, however, a true distinction be- 
tween an actual critical (global) temperature T c and a or- 
dering temperature T* is difficult to achieve. This is not 
unlike the problem encountered in the strong-coupling 
limit, where, in the context of superconductivity, Cooper 
pairs form at some large T* , but condense in a Kosterlitz- 
Thoulcss transition only at a T c much lower. Neverthe- 
less, the ideas presented here have more in common with 
those discussed, e.g., in diluted magnetic semiconductors, 
where, in the limit of small impurity doping, one may ob- 
serve the formation of isolated, locally ordered clusters 
at a large temperature T*; as the temperature is lowered 
the associated wavefunctions expand and start to over- 
lap at T c in a percolative transition^. We assume that 
this scenario gives the best description of the SC tran- 
sition in the undcrdoped regime. Therefore, we adopt 
the viewpoint that a phase transition in the inhomoge- 
neous regime occurs once a certain percentage of lattice 
sites assumes a finite-value of the OP, either AF or SC, 
and chose here as the critical percentage p c =60%. This 
value is motivated by extensive MC calculations for site- 
diluted 2D Hciscnbcrg models that have shown that the 
percolation transition (at T=0) is essentially a classical 
one, i.e., long-range order (LRO) is lost once ~40% of the 
spin-carrying sites are removed, in agreement with (clas- 
sical) bond-diluted random resistor networksS*2&. Since 
the SC phase transitions are in the same (XY) univer- 
sality class, we adopt this particular value for that case 
as well. Note that because of the previously discussed 
proximity effects the relative amount of sites carrying a 
finite OP is actually larger than the number of sites with 
a finite value of U or V in Huf', this is the main effect of 
the hopping terms in the Hamiltonian when compared to 
the localized spin models employed in MC simulations. 
Even if the actual p c differs from the chosen value - and 
this can only be verified by a more thorough investiga- 
tion into site-or bond-diluted superconductors, presum- 
ably with more refined MC simulations - the shape of 
the T c (x) curves is not supposed to be extraordinarily 
affected by any particular choice of p c . 

Unfortunately, it turns out that even for a strongly 
doped system the AF OP in isolated pockets survives al- 
most up to T^ F , probably an unfortunate consequence 
of the approximation involved. Therefore, we have in 
addition calculated the approximate temperature upon 
where the associated gap in N(ui) closes and defined this 
value as T*. This procedure should not be regarded as 
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FIG. 11: Temperature-dependent N(u>) for two different dop- 
ing levels (n)=0.90 (a) and (n)=0.80 (b). Whereas in (a) the 
gap-like feature survives almost up to the clean transition 
temperature T c F ~l, it is barely notifiable in (b) at temper- 
atures just above T c . 



particularly unusual, since (a) it is a common way to 
determine T* in underdoped cuprates anyway, and (b) 
because T* in itself is not a strictly defined quantity. 
Rather, it is a temperature much larger than T c below 
which certain excitations are suppressed, and the actual 
values cited in the literature vary depending on the ex- 
perimental quantity investigated, or even allow for the 
definition of multiple PG temperatures. 

The resulting phase diagram (for the set of parameters 
employed throughout this work) is shown in Fig ll2l It 
is characterized by a multitude of different relevant tem- 
peratures in the underdoped regime, such as the critical 
(SC) temperature T c (thick broken line), a temperature T 
up to where local superconductivity survives (indicated 
by stars in Fig ll2f> . and a PG temperature T*, deter- 
mined by the gap in N(u). Most striking is the sim- 
ilarity to the standard cuprate phase diagram, with a 
PG phase and the SC dome and it is noteworthy that 
this has been achieved within such a simple scheme and 
without any further assumptions. In the framework of 
the mixed-phase state as discussed here, the central role 
is assumed by the level of optimal doping, which sepa- 
rates the inhomogeneous and homogeneous parts of the 
phase diagram from each other. For doping levels smaller 
than x c , doping is tantamount to replacing AF-ordered 
sites with SC ones, whereas the density within those SC 
clusters remains unchanged (at least in a first approxi- 
mation). For x>x c , on the other hand, further doping 
actually does lead to a decrease of the charge density 
across the whole sample, and subsequently to a reduction 
in T c . The existence of a maximum in T c is therefore a 
natural consequence in such a scenario. It needs to be 
noted, however, that according to Fig^|the reduction in 
T c is much smaller in the overdoped state than what is 
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FIG. 12: Phase diagram for the cuprates, derived from 
Eq.Q, where hole doping leads to a coexistence of SC clus- 
ters (density nsc^l-x c ) with the parent Mott insulator. The 
line of optimal doping x c separates two fundamentally dif- 
ferent regimes, the homogeneous (h) overdoped phase, where 
rather common BCS physics can be expected, and the in- 
homogenous (ih) underdoped regime, which consists of the 
AF/SC mixture and whose SC transition is of percolative na- 
ture. T c determined as described in the text, lines serve as 
guide to the eyes. Local SC clusters survive up to T, which 
forms the upper bound of a SC fluctuating (SC/F) regime, 
and T* denotes the PG temperature as derived from the den- 
sity of states. The existence of the comparatively large gap for 
small doping levels allows for the definition of an additional 
characteristic temperature, if so desired. This temperature, 
however, is poorly defined and thus omitted here. We have 
also indicated the boundaries of the AF phase. Assuming that 
V is in some sense derived from spin fluctuations, the T c (n) 
curve in the overdoped regime is replaced by the thin bro- 
ken line, and similarities to the conventional phase diagram 
become more pronounced. 



actually observed. However, this calculation rests on the 
assumption that the interaction V does not depend on x, 
whereas in reality (if, for example, AF spin fluctuations 
are considered) V presumably has to be reduced in par- 
allel to (n). We have indicated this by the thin- broken 
line in Fig ll2l which shows T c assuming V r oc(n/nsc) 2 , 
as one would expect for an interaction mediated by spin 
fluctuations. 

We also want to mention that T is surprisingly large 
even for £=0.05, where it is about half the value of T c . 
Detecting T may not be so easy, however. The excita- 
tion spectrum does not change appreciably upon crossing 
T - the FS along the BZ, already established, remains 
unaffected, whereas along the anti-nodal direction the 
proximity-induced gap in all likelihood suffices to prevent 
the creation of a FS, which would otherwise occur. At 
best, one can expect a small increase in spectral weight at 
E-p as this temperature barrier is crossed. The tempera- 
ture regime between T and T c (SC/F) can with justifica- 
tion be called a "fluctuating" regime. This is terminology 
borrowed from recent investigations into the Ncrnst effect 
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in LSCO and YBCG2LS, leading to the establishment of 
an additional temperature T u , below which vortex-like 
excitations can be observed, and thus local (preformed) 
superconductivity may exist. Although T v has the same 
doping dependence as T (as is found here), it reaches 
more than 100K in optimally doped LSCO, in apparent 
disagreement with the phase diagram above. Neverthe- 
less, the question of an exact definition of T v is open as 
is a full investigation of the effects of magnetic fields on 
mixed AF/SC regimes still lacking. The phase diagram of 
Fig ll2l is also in remarkable agreement with one proposed 
recently for La-doped Bi-2212Si. After a careful analy- 
sis of resistance data, there the authors demonstrated 
that in the regime immediately above T c , the resistance 
can be described via a simple universal two-component 
scaling function, one component with an insulator-like 
weak temperature dependence and comparably large re- 
sistance (i.e., the hole-poor phase), and a second metallic 
one with small resistance and related (by the authors) 
to SC fluctuations. This is precisely the scenario de- 
scribed here. Resistor network calculations similar to 
those sometimes used in manganites may help to confirm 
these results for the model considered here. NMR in- 
vestigations for underdoped LSCO^ paint a very similar 
picture by demonstrating that the hole density x\ varies 
on a nm length scale. The reported standard deviations 
in the local charge density, Axi, are comparable with the 
ones found here if both the temperature and the assump- 
tions of model are taken into account. Finally, we also 
want to mention here that the picture emerging for the 
SC part of the underdoped phase is closely related to the 
phase-separation scenario proposed by Emery and Kivcl- 
son more than a decade ago^i. Furthermore, the exis- 
tence of more than one characteristic temperature above 
T c in the underdoped regime was recently discussed by 
Lee et alw^. 

Naturally, Tg C =T c in the overdoped limit, since it is 
regarded as a homogeneous system, governed by con- 
ventional BCS physics. Interestingly, the fact that p c 
is larger than 50% immediately suggests the opening of a 
"glassy" phase without LRO, which will occur once AF 
and SC volume fractions are approximately balanced. If 
quantum fluctuations were involved, they would presum- 
ably tend to increase the percolation threshold, i.e., sta- 
bilize the phase without LRO. Yet, this should only be 
a minor effect and leave the generic phase diagram un- 
changed. Nevertheless, due to the discussed proximity 
effects, an SC/AF coexisting phase may appear as well 
and is actually present and indicated in FiglT2l 

If one assumes that doping in cuprates in general is fa- 
cilitated by introducing hole-rich clusters with (in a first 
approximation) x c , than xaf will be proportional to x c . 
In fact, the fairly quick disappearance of the AF phase 
is a logical consequence, for it is caused by the compar- 
atively large size of the metallic/SC clusters, and which 
invalidates the naive expectation of an AF percolation 
threshold xaf~30%-50%^. The key issue here is that 
it is the partial derealization of holes that is responsible 



for the much quicker than anticipated loss of antifcrro- 
magnetism. For example, for the realistic value x c ~0.15 
(rather than x c used here) and assuming that antifer- 
romagnetism vanishes according to the rules of classical 
percolation, £af=0.06, a value close to the experimen- 
tally realized one. In parallel, the onset of SC order is 
predicted for xsc=0-09, with a glassy regime in between. 
Nevertheless, these considerations only hold for a truly 
2D system; if, for example, a bilayer material such as Bi- 
2212 is considered, the transition into the SC state will 
take place at smaller doping values since the percolation 
threshold is reduced in 2D when compared to layered or 
even truly 3D systems, assuming here that the value of x c 
is not affected by the lattice geometry. This consideration 
also applies to the AF phase, and thus bi- and tri-layered 
materials should almost always prefer a SC/AF coexist- 
ing regime rather than the spin-glass phase encountered 
in single-layer materials. 

Finally, we wish to comment on another subtle issue. 
The reader may have noticed that our implementation 
of the doping process is not entirely consistent. If dop- 
ing is really a random process, and there are no indica- 
tions otherwise, the charge-depleted volume fraction will 
not increase linearily with doping. Rather, clusters will 
sooner or later start to overlap, and thus the increase 
in metallic volume will be slower than linear. Conse- 
quently, the density of holes in those areas will increase 
somewhat from the original level, here assumed as the 
one of optimal doping, because the holes are squeezed 
into a relatively smaller volume. Therfore, on should as- 
sume that initially the hole-density of metallic neighbor- 
hoods is smaller than x c . If, for example, originally the 
charge density of metallic clusters were ngc^O.90, then 
this would in parallel shift the breakdown of antifcrro- 
magnctism to xaf=0.04 rather then the value of 0.06 
introduced above^i. x c would still be the doping level 
upon where the system becomes homogeneous, which has 
to happen eventually, and the relation between xaf and 
x c could be figured out by geometric considerations, for 
example by covering the lattice with randomly placed 
circles of a certain radius, and determining when the un- 
derlying lattice is fully covered^. The decrease of density 
in hole-rich cluster is then reflected in the observed small 
downward-shift of the chemical potential within the MH 
gap, which was remarked upon in Section III and ob- 
served experimentally^. 



VI. CONCLUSIONS 

Motivated by a series of recent PE experiments we have 
explored the physics of underdoped cuprates under the 
assumption that they represent a mixture of the undoped 
parent insulator and a SC/metallic phase with a hole- 
density similar to the one of optimal doping. As empha- 
sized by us previously^, this scenario is consistent with 
most of the PE data, specifically the two-branched signals 
that are observed throughout the underdoped regime, as 
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well as the development of the nodal FS arc with doping. 
Disorder as it manifests itself in spatially varying chemi- 
cal potentials seems to be the driving force in LSCO, but 
inhomogeneity is not necessarily tied to disorder alone. 
Here, we point out the importance of proximity effects for 
such exotic electronic states, and how they can explain 
even subtle details of experimental results, such as the 
development of the low-energy excitation gap in the anti- 
nodal direction and the behavior of the spectral weight 
at Ep. Proximity effects are always present in inhomoge- 
neous electronic systems, although they may be hidden if 
both competing orders are of sufficiently equal strength. 
We have also extended our previous calculations to fi- 
nite temperatures and predict a gradual increase of the 
FS arc for temperatures above T c and that the intensity 
along the BZ diagonal should remain constant, both in 
principle verifiable directly by ARPES measurements. 

In the underdoped phase, doping is tantamount to in- 
serting clusters of hole density x<x c into the insulating 
host. The accompanying phase transitions are percola- 
tive in nature, irrespective of whether they are caused by 
a change in temperature or doping. Pseudogaps appear 
as hallmarks of remnant antifcrromagnctism and mostly 
arise because of the vast difference in energy between 
SC and AF excitations. Transport properties aside, the 
superconductor-metal transition for x<x c is a smooth 
one, quite unlike what is observed for overdoped mate- 
rials, where it is marked by the sudden appearance of 
Bogoliubov quasiparticles in both experiment and the- 
ory. Above T c , the underdoped phase possesses isolated 
SC clusters, which may survive even into the AF phase 
and up to surprisingly large temperatures. This gives 
rise to an additional characteristic temperature T, as in- 
dependently suggested by Nernst effect and resistance 
measurements. In addition, the variation in local hole 
density is approximately as large as the one recently re- 
ported in LSCO. 

Overall, it was shown that those fairly simple consider- 
ations can explain a substantial part of the high-T c phe- 
nomenology, including the phase diagram. Moreover, no 
further assumptions of competing states, changes in pa- 
rameters, curious non-Fermi-liquid metallic states, etc. 
have to be made to rationalize an enormous variety of 
experimental data. From this point of view, the most im- 
portant task in cuprate theory is to decipher the effects 
of chemical doping, particularly with respect to possible 
spatial chemical potential fluctuations. 

We also want to stress here that similar ideas 



have been widely discussed in manganites^, and re- 
cent ARPES data on the doped bilayer material 
La2_2a;Sri +a ;Mn2 07 are remarkably similar to those of 
underdoped cuprates— Most importantly, a two-peak 
structure in spectral intensity was found along the BZ 
diagonal (though not along the anti-nodal direction) for 
x=0A. Since we have shown that such features are in- 
evitable consequences of mixed-phase states, it seems 
natural to interpret those data along the same lines, the 
only difference being the competing phases in question. 
In manganites, one of the phases is the ferromagnetic 
metal, whereas the competing insulating one is either 
a CO one and/or an AF insulator. Nevertheless, given 
what we know about the cuprates, this explanation seems 
more reasonable and natural than the one suggested by 
the authors of RefliSl which relied on strong electron- 
phonon coupling, especially since the evidence for a phase 
mixture is stronger in manganites than in cuprates to 
begin with—. The underlying theme of emerging com- 
plexity in both cup rates and manganites was recently 
discussed in Refl49t 

In this context, it would be particularly interesting to 
obtain ARPES results for slightly doped Lai-zSr^MnOa, 
since this material is most closely related to LSCO and 
has a comparatively simple phase diagram, with only the 
half-filled Mott insulator and the ferromagnetic metal- 
lic phase present for doping levels smaller than x=0.5. 
If the mixed-state picture applies there as well, the re- 
sulting PE spectra will look just like the ones in LSCO, 
including a PG as postulated earlier—. In particular, 
based on the proximity effects emphasized in this work, 
one would expect an excitation gap in the anti-nodal di- 
rection even though the ferromagnetic metallic phase is 
nominally ungapped, unlike the SC one. From what we 
know about spectral functions in electronically inhomo- 
geneous systems, such experimental results, if they be- 
came available, should be regarded as the best indicator 
of phase separation in slightly doped manganites, ruling 
out the existence of a canted phase. 
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